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Abstract 

One-dimensional frost growth and densification in laminar flow over flat surfaces has 
been theoretically investigated. Improved representations of frost density and effective 
thermal conductivity applicable to a wide range of frost circumstances have been 
incorporated. The validity of the proposed model considering heat and mass diffusion in the 
frost layer is tested by a comparison of the predictions with data from various investigators 
for frost parameters including frost thickness, frost surface temperature, frost density and 
heat flux. The test conditions cover a range of wall temperature, air humidity ratio, air 
velocity, and air temperature, and the effect of these variables on the frost parameters has 
been exemplified. Satisfactory agreement is achieved between the model predictions and the 
various test data considered. The prevailing uncertainties concerning the role air velocity and 
air temperature on frost development have been elucidated. It is concluded that that for flat 
surfaces increases in air velocity have no appreciable effect on frost thickness but contribute 
to significant frost densification, while increase in air temperatures results in a slight increase 
the frost thickness and appreciable frost densification. 
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Nomenclature 
b plate thickness (m) 
c p specific heat (kJ/kg K) 

D mass diffusion coefficient (m 2 /s) 
h c convective heat transfer coefficient (W/m 2 K) 

h m convective mass transfer coefficient (kg/m 2 s) 
k thermal conductivity (W/m K) 

k„ dc effective thermal conductivity of air with mass diffusion and eddy convection (W/m K) 
k f effective thermal conductivity of frost (W/m K) 

L plate length (m) 

L sv latent heat of sublimation (kJ/kg) 

M molecular weight (kg/mol) 

Nu Nusselt number, h c Llk 
p pressure (Pa) 

p v partial pressure of water vapor at temperature T (p a ) 
p v0 partial pressure of water vapor at temperature Tq (Pa) 

Pr Prandtl number 
q heat flux (W/ m 2 ) 

R universal gas constant (kJ/kmol K) 

R a gas constant for air (kJ/kg K) 

R v gas constant for water vapor (kJ/kg K) 

T temperature (K) 

T f frost temperature (K) 

T„ freestream temperature (K) 
u m freestream velocity (m/s) 

W plate width (m) 

x coordinate normal to the plate, i.e. along the frost thickness (m) 
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Greek Symbols 

<j> relative humidity (%) 

^ porosity (void fraction) 

v kinematic viscosity (m 1 2 /s) 

p j- frost density (kg/m 3 ) 

C k mtc / k jcc 

co humidity ratio of moist air, kg (water vapor)/kg(dry air) or kg v /kg a 

X solids (ice) fraction 

Subscripts 

0 reference 

1 continuous medium (stagnant gas) 

2 solid particle (dispersed) 

oo ambient (ffeestream) 

a dry air 

c eddy convection, critical 

d diffusion 

e effective medium 

/ frost 

m melting 

s frost surface 

t total (ambient) 

v water vapor 

1. Introduction 

Frost formation represents an important consideration in cryogenics, refrigeration, 
aerospace, meteorology and various process industries. Processes involving simultaneous 
heat transfer and frost formation are occasioned in gas coolers, refrigerators, regenerators, 
freeze-out purification of gases, cryopumping, and the storage of cryogenic liquids. 

In reality, frost offers significant thermal resistance to heat transport in heat exchangers, 
and increases the pressure drop in channels by flow restriction. Frost/icing on aircraft wings 
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causes significant aerodynamic penalties in lift and drag, leading to a decrease in lift and an 
increase in drag. Frost/ice formation on cold fuel/oxidizer tanks and other components may 
cause debris concerns when it is shed on critical components of launch vehicles. These 
applications call for an accurate prediction of frost growth, frost density, frost surface 
temperature, and total heat flux. A fundamental understanding of the nature of frost 
formation and simultaneous heat and mass transfer during frost growth is thus of great 
practical importance. 

Frost forms when moist (humid) air comes in contact with a cold surface whose 
temperature is less than the freezing temperature of water (273 K), and is less than the dew 
point temperature, so that water vapor directly passes from a gaseous state to a solid state. A 
schematic of the frost layer is shown in Fig. 1 . During frost formation, a portion of the total 
mass flux of water vapor diffuses through the frost layer, and undergoes phase transformation 
to ice which deposits on the existing frost crystals, resulting in frost densification (Jones and 
Parker [1]; Na and Webb [2]). The remainder of vapor mass flux freezes at the frost surface, 
and serves to increase the frost thickness. A review of frost formation in simple geometries 
(flat plates, parallel plates, cylinders, and annuli) was provided by O'Neal and Tree [3]. 

Experimental studies for frost formation on flat surfaces were published by various 
investigators; see Jones and Parker [1], Yonko and Sepsy [4], Ostin and Anderson [5], Lee et 
al. [6,7], Cheng and Wu [8], Hermes et al. [9] among others. These studies suggest that 
during frost formation, frost thickness, frost surface temperature, frost density, and frost 
thermal conductivity all increase continuously. Both the density and the thermal 
conductivity are found to increase by an order of magnitude. As the frost layer thickens, the 
insulating effect of the frost layer predominates resulting in an increased thermal resistance 
of the frost layer and a decrease in heat transfer rate (Sherif et al. [10]). Measurements also 
reveal that the frost growth, densification and heat flux depend on a number of parameters 
such as frost density, frost surface temperature, plate surface temperature, and environmental 
parameters (air humidity co^, air velocity u^ y air temperature , etc.). Experiments by 

Ostin and Anderson [5] and others suggest that the frost growth depends strongly on the cold 
plate temperature and the air humidity. 

A number of theoretical models for frost formation on flat surfaces have been proposed 
for frost formation on flat surfaces. Notable among these models include those of Brian et al. 
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[11,12], Jones and Parker [1], Sami and Duong [13], Sheri fet al. [10], Tao et al. [14], Sahin 
[15], Le Gall et al. [16], Lee et al. [6], Cheng and Cheng [17], Na and Webb [2], Hermes et 
al. [9J. These models differ primarily in the treatment of frost density and thermal 
conductivity correlations and/or models and of mass diffusion within the frost layer. A 
comprehensive review of various theoretical models is provided by Tao et al. [14]. 

Notwithstanding the appearance of a large number of articles on frost growth 
measurements and modeling, there exist some inconsistencies and uncertainties primarily 
with regard to the role of air velocity and air temperature on frost growth and densification. 
Yonko and Sepsy [4] report that for flat surfaces even though increasing air velocity has no 
appreciable effect on frost thickness, the frost density is considerably increased. Data by Lee 
et al. [6] for flat plates suggest that increasing air velocity slightly increases the frost 
thickness, but changes the frost surface temperature (and thus frost density) appreciably. 
Measurements by Shah [18] on a circular plate indicate that higher Reynolds numbers and/or 
higher gas temperatures yield higher frost density and a decrease in frost thickness (Sami and 
Duong [13]). Data by Cheng and Wu [8] for flat plates indicate that increases in air velocity 
or air temperature resulted in little change in frost thickness but produced a denser frost layer 
on the basis of photographic (visual) examination. 

In view of the aforementioned uncertainties in frost formation, it is evident that a more 
detailed understanding of the effects of frost parameters on the frost growth and 
densification is requisite. It is the purpose of this work to further clarify the role of air 
velocity and air temperature on the frost growth and densification processes. This objective is 
carried out with the aid of a one-dimensional theoretical model proposed to predict the frost 
growth and densification in laminar flow over flat surfaces. Improved formulations for frost 
density and frost conductivity valid for a wide range of conditions are incorporated. The 
validity of the proposed model will be assessed by comparison with measurements frost 
parameters includingyra$7 growth, frost surface temperature, frost density and frost heat flux 
as reported by various investigators covering a range of test conditions including cold plate 
temperature , air humidity , air velocity and air temperature. 
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2. Proposed Model 

In the present work, we propose a one-dimensional model for frost growth and 
densification along the lines of Jones and Parker [1 ] and Cheng and Cheng [17], Improved 
formulations of frost density and thermal conductivity, applicable for a wide range of frost 
circumstances, are incorporated, and this marks an important difference between the present 
model and those of other investigators. 

2. 1. Physical Assumptions 

The following physical assumptions are made in the present analysis: 

(i) Heat and mass diffusion within the frost layer are assumed quasi-steady and one- 
dimensional. 

(ii) The frost layer grows in a direction normal to the cold plate surface with uniform frost 
thickness along the plate length. 

(iii) Frost density and thermal conductivity are independent of frost thickness, and represent 
values averaged over the frost layer thickness. Thus at any instant, they are thus assumed 
uniform within the frost layer. Measurements by Brian et al. [12] suggest that at any instant 
during frost growth density variations within the frost layer are not significant. 

(iv) The air pressure is uniform in the airstream and within the frost layer. Hermes et al. [9]. 

(v) The Lewis analogy for heat and mass transfer is applicable, and Lewis number of unity is 
considered. Measurements and analysis by Lee et al. [6] indicate that Lewis number 
increases with air velocity and air humidity, and is in the range of 0.7 to 1.0. 

(vi) Both the heat and mass transfer coefficients at the frost surface are considered constant. 

(vii) The gas spaces in the frost layer are small enough that heat transfer by natural 
convection within the frost layer is negligible (Woodside [19]). Thermal radiation within the 
frost layer is negligible, considering that the frost temperatures are low (Dietenberger [20]). 
(viii) The total gas pressure is constant within the frost layer, and is equal to the atmospheric 
pressure (Le Gall et al. [16)). 

(ix) Soret and Dufour effects are negligible (Le Gall et al. [16]). 

(x) The local humidity of air in the frost layer corresponds to saturation humidity at that 
temperature (thermodynamic equilibrium). 
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(xi) It is also assumed that at the frost surface the water vapor is saturated. Strictly speaking 
this assumption is somewhat tenuous, although many frost growth models consider this 
assumption. Detailed investigations by Na and Webb [2] indicate that water vapor is 
generally supersaturated at the frost surface rather than saturated, leading to appreciable 
overprediction of the mass transfer rate of water vapor onto the frost surface. The 
supersaturation degree (measure of departure of the actual water vapor pressure at the frost 
surface to the saturated water vapor pressure at the frost surface) is strongly dependent on the 
surface energy, which is a function of the surface coating governing the contact angle. In 
view of the absence of such information in reported experimental data, we follow for the 
present purposes the simplifying assumption of saturated vapor condition at the surface. 

2.2. Frost Density and Effective Thermal Conductivity 

The prediction of frost growth requires information concerning the average frost density 
and effective frost thermal conductivity, both of which depend on the porosity of the frost 
layer. The frost porosity y is directly related to the frost density pj- by the relation 

^ = (l-P//picc)/(l-P g /Acc) (la) 

where /? ice is the density of ice (usually taken as 917 kg/m 3 ), and p g the density of air- vapor 
mixture. The porosity is related to the solids (ice) fraction x by the relation 
X = \-'F (lb) 

A detailed review of frost density and thermal conductivity correlations (and their range 
of applicability) is provided by Tao et al. [14], O'Neal and Tree [3], Shin et al. [21], and 
Yang and Lee [22] presented a comprehensive list of the correlations for the effective 
thermal conductivity of frost. 

2.2.1. Frost Density 

Recently, the author [23] reviewed the existing data and correlations for frost density, and 
proposed a new and improved correlation in terms of dimensionless frost surface temperature 
and flow Reynolds number for application to laminar forced flow over flat surfaces. This 
correlation is based on simultaneous measurements of frost density and frost surface 
temperature and is shown to be valid over a broad range of frost density, and is expressed by 
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(2a) 


— = O.50exp{- [0.376 + 1 ,5(l - 0)](l - £° 5 )} 

/^ice 

where 

° = (T fi -T w W m -T w \ £ = Re/ Re r (2b) 

In the preceding equation, d? is the dimensionless frost surface temperature, and Re c is 

the critical Reynolds number for laminar-turbulent transition. The Reynolds number is 
defined by 

Re = u a L / v ( 2 C ) 

where L is the characteristic plate length, w M is the air velocity, and v m is the kinematic 
viscosity of the air. Fig. 2 (adapted from [23]), displays the variation of frost porosity and 
density as a function of dimensionless frost surface temperature with Reynolds number as a 
parameter. This model was validated over the full range of frost surface temperature. Note 
that the widely used correlation of Hayashi et al. [24] expresses frost density dependent on 
frost surface temperature only, and is known to considerably overpredict the frost density at 
relatively high frost surface temperatures close to the melting temperature [14]. 


2.2.2. Effective Frost Thermal Conductivity 

The effective thermal conductivity of frost k f is obtained from the work of the author 

[25] on the basis of a porous packed bed model of Zehner-Schlunder [26] considering a unit 
cell approach and one-dimensional heat flow (see Hsu etal. [27], Cheng and Hsu [28]): 
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where 


(3a) 


C = k uJc /k ice (3b) 

In the preceding equations, k mle denotes the effective (stagnant) thermal conductivity of 

entrapped air considering mass diffusion and eddy convection. The shape factor B is related 
to the porosity (defined by Eq. 1) by the relation 

B ~ C[{\-y)ly} m (3 C ) 
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with in = 10/9 . The constant C is taken as 2.5 for cylindrical particles (Zehner-Schlunder 
[26]), which were found to represent frost structure and resulting frost thermal conductivity 
satisfactorily [25], A detailed comparison of Zehner-Schlunder model with data (for spherical 
particles) for ice to air thermal conductivity ratio ranging from 8 to 1200 and for various 
values of solids fraction was recently presented by the author [29]. 

In the presence of mass diffusion and eddy convection, the effective thermal conductivity 
of air k mk is expressed as (Kandula [25]) 

k «* = k « + k „ + K (3d) 

where k u is the molecular thermal conductivity of air, kj is the contribution of mass diffusion 


to the effective thermal conductivity of air, and £ c denotes the contribution of eddy 
convection to the effective conductivity of air. 

The mass diffusion contribution k d Eq. (4) is obtained as [25] 


= L 2 d{mjm v ) 


R]T} 


PvO eX P 



f---T 

K 

Jo T f)_ 


(3e) 


In the preceding equation, D is the mass diffusivity of water vapor, L sv is the latent heat of 


sublimation, M a ,M y denote the molecular weights of air and water vapor respectively, p v0 
is the partial pressure of water vapor at reference temperature T 0 , R v is the gas constant for 


water vapor, and T f is the frost temperature. 

The effect of eddy convection [30] prevalent in the air pockets is shown to be 
satisfactorily represented by the condition [25] 

k c lk m *\ (30 

From a more fundamental consideration the effects of eddy convection are intimately 
related to the phenomenon of thermal dispersion (see Cheng and Vortmeyer [31]; Hsu and 
Cheng [32]). Convection in variable porosity media, including channeling effects, were 
considered by Vafai [33]. 

Thus the effective thermal conductivity of frost is obtained from Eqs. (3a-30- Fig. 3 
(adapted from [25]) shows a comparison of the predicted thermal conductivity with the well- 
established correlations of Sturm et al. [34], Van Dusen [35], and Ostin and Anderson [5], 
More detailed comparisons of the predicted thermal conductivity with the measurements of 
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Pitman and Zuckerman [36], Brian et al. [36] and others at various temperatures are 
presented in the author's work [25], 


2.3. Governing Equations and Boundary Conditions 
2.3. L Governing Equations 

An energy balance for an elemental frost thickness in the interior of the frost layer 
(considering latent heat transport from the water vapor) can be expressed as [1], see Fig. 1): 


Al 


dm „ 


dx 1 JV dx 

where m x stands for the local mass flux of water vapor in the x direction (normal to the 


(4a) 


plate). The quantity m x is given by the mass diffusion equation 


m x = -D 


d P, 

dx 


(4b) 


where p v stands for the local density of water vapor in the frost layer, so that 

dlh * _ D ^P V 
dx dx 2 


(4c) 


2. 3. 2. Boundary and Initial Conditions 
a. Mass Flux at the Frost Surface 

The total mass flux of water vapor transferred from moist air to the frost surface can be 
written as 

"h = -a.) (5a) 

where h m denotes the convective mass transfer coefficient. The quantities co u , co s are 

respectively the humidity ratios in the atmospheric air and at the frost surface. Eq. (5a) can be 
recast as (Jones and Parker [1]) 


d t \ dx s dPf 

a**. dt 


(5b) 

H 

II 

h 

CD 

(5c) 


Eq. (5b) physically implies that a portion of the total mass flux of water vapor is diffused into 
the existing frost layer before it freezes, and serves to increase the frost density. The 
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remainder of this vapor is transported to the frost surface and freezes at the frost surface, and 
serves to increases the frost layer thickness (Jones and Parker [1 ]; Na and Webb [2]). 

b. Heat Flux at the Frost Surface 

The total heat(energy) transferred from the air to the frost surface is comprised of two 
components: sensible heat transfer (by convection) by virtue of the temperature difference, 
and the latent heat of sublimation due to humidity difference. Mathematically, 

<7, = <7 sen + 4 In, (6a) 

where 

<7 sen = K i r „ ~T fi \ <7 la, = '»,4v (6b) 


In the preceding equation, m t is given by Eq. (5b), and h c is the convective heat transfer 
coefficient. Eq. (6a) can be recast as 




4 - X 


d£j 

dt 


(6c) 


The quantity in the parenthesis represents the energy crossing the frost surface by conduction 
(Cheng and Cheng [17), so that 


9cond ^1 


( dT} 
\dx J s 


dx 


M 


c. Temperature Boundary Conditions at the Wall and the Frost Surface 
The boundary temperatures are prescribed as 
T = T. x = 0 


T = T, 


x = x. 


(7a) 

(7b) 


d. Initial Condition 

The initial condition is specified by the requirement that the frost thickness is zero. Thus 
x s — 0 t = 0 (8) 


2. 4. Temperature Distribution in the Frost Layer 

By virtue of the assumption that the frost density remains uniform across the frost layer 
thickness, the amount of water vapor being frozen in the interior must be the same at all 
locations in the frost layer (Jones and Parker [1]). Consequently, Eq. (4c) becomes 
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dm x _ m d 
dx x s 

Thus Eq. (4a) governing the temperature field is expressed by 
. d 2 T / - 


( 9 ) 


f dx 2 


= -L 


m, 


\ x > 


(14) 


Substituting for m d from Eq. (5c), and integrating twice with the appropriate boundary 

conditions (Eqs. 6d, 7a and 7b), we finally arrive at the temperature distribution in the frost 
layer as follows: 

T = T w + c { x + c 2 x 2 


K f 


c,=- 


k f X s 


K\®a-a s )-p f ^- 


(11a) 

(lib) 

(11c) 


The temperature at the frost surface corresponding to x = x s can at once be obtained from 
Eq. (1 la) as follows (Cheng and Cheng [17]): 


T = T.+ 


'x ' 


\ k f J 




v k f ) 


-aJ+Ps ~ 


( 12 ) 


2. 5. Frost Layer Growth Rate 

From Eqs. (5a) and (5b), the frost layer growth rate is obtained as 
. ( \ dx dp. 

On the basis of the correlation for the frost density (Eq. 2), we have 
Pf 1 P icc = O.50exp{[c, + c 2 (l - 0)](l - £ 05 )} 


Thus 

dpj_ _ dp f d9 dT fs dx 2 
dt dO dT fs dx s dt 

Performing the differentiation and rearranging, we finally obtain that 


(13a) 


(13b) 


(13c) 
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(13d) 


(T m ~T w )dx, J 


(14a) 


2. 6. Heat and Mass Transfer Coefficients 

The convective heat transfer coefficient h c for laminar force flow over a flat plate is 
determined from the well-known correlation (Kays and Crawford [37]) 

Nu - 0.664 Re l/2 Pr l/3 

where the Nusselt number Nu and the Prandtl number Pr are defined by 

Nu = h c L/k, Pr = c p /atk (14b) 

and Prefers to the thermal conductivity of air The fluid properties for the boundary layer are 
evaluated at the average film temperature. 

The convective mass transfer coefficient h m is obtained on the basis of Chilton-Colbum 

analogy between heat and mass transfer (Eckert and Drake [38]), which is found to be 
applicable under frosting conditions (O'Neal and Tree [5]): 

ShLe 2n = Nu (15a) 

where the Sherwood number Sh and the Lewis number Le are defined by 
Sh = h m LI D, Le = a / D = Pr/ Sc (15b) 

In the preceding equation, a - k/ pc p represents the thermal diffusivity. Sc = Dl v is the 
Schmidt number, and c p is the specific heat of moist air. Consequently we have 
h n,= h c l (p c p LeVi ) (15c) 


with specific heat of moist air quantity c determined from the mixture relation 


c - c + CO c 

p pa a pv 


(15d) 


A value of Le = 1 is considered in the present analysis. 
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2. 7. Thermophysical Properties 

The thermal conductivity of air is calculated as (Rohsenovv and Hartnett [42]; 
Dietenberger [20]) 


k... 


2.646x1 0"' 3 * 


’ 1/2 


[l + (245/r)l0 


12/7 


w 

mJC 


( 16 ) 


where T is in degrees Kelvin. The thermal conductivity of ice is obtained from the relation 
(Dillard and Timmerhaus [43]; also Dietenberger [20]) 

k icc =630/7’ W/(m.K) (17) 

where T is the absolute temperature (degrees K). 

The mass diffusivity is determined from (Pruppachar and Eklett [39]) as 
D = 2.1 lxl 0" 5 (7^ /r o )' 94 m 2 /s (18) 


where Tp is the frost temperature in degrees Kelvin, and Tq the reference temperature 
(273.15 K). 

The latent heat of sublimation of ice is determined from (Mago et al. [40]) 

= {-0. 1 083[l . 8(r - 273. 1 6) + 32]+ 2833} kJ/kg (19) 

where T is in degrees Kelvin. 

The humidity ratio of air at the frost surface &> s is determined from (Kays and Crawford 
[37]): 

co s = 0. 622/7 vJ /(p t - p vs ) (20) 

where p t is the total pressure in atm, p vs is the saturation vapor pressure corresponding to the 
frost surface temperature T fs . The vapor pressure p vs is evaluated from standard curve-fits 
(Mago et al. [40]; Flatau [41]). 


3. Numerical Solution 

The frost growth and frost surface temperature are obtained by numerically integrating 
Eq. (13d) in conjunction with Eqs. (6), (12) and (13b), considering the initial condition 

expressed by Eq. (8 ). The equation for the frost surface temperature is implicit. Initially a 
frost surface temperature at t + At is guessed, and a bisection method is used to iterate for the 
frost surface temperature. A time step of A/ = 0. 1 sec is chosen. For stability of computation, 
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a very small value of x s =10 5 m is considered. The calculation is terminated when the frost 

surface reaches the melting temperature. From the calculated first surface temperature, the 
frost density is evaluated with the aid of Eq. (2). The heat flux values (sensible, latent and 
total) are subsequently obtained from Eq. (6a) and (6b). 

The numerical solution is found to be robust over a wide range of test parameters 
considered in the present investigation. 

4. Results and Comparison 

The validity of the proposed model is tested by a comparison of the predictions of frost 
properties with test data from various investigators comprising a range of frost conditions 
(cold plate temperature, ambient air humidity, air velocity, and air temperature). All the test 
data chosen for comparison correspond to the geometry of flat surfaces. Although frost 
thickness is measured in every test data, only certain data sets contain frost surface 
temperature, frost density, and heat flux. Note that accurate frost surface temperature 
measurement is relatively difficult. In all the comparisons, a value of Re c % 10 5 appearing in 
the frost density correlation, Eq. (2), has been considered for the model predictions [23]. 

4. 1. Comparison with Data ofYonko and Sepsy [4] 

Figure 4a-4c displays a comparison of model predictions for frost thickness with the 
experimental data ofYonko and Sepsy [4] on an aluminum plate. The test plate dimensions 
are L =0.3398 m, W= 0.1905 m, and b =0.0222 m. The frost thickness is measured to 
within 0.13 mm with the aid of a micrometer working against a spring-loaded cross-slide. 
Uncertainly in the measured frost thickness is about 0.13 mm. These data represent one of 
the early measurements on frost formation, and have been very widely considered in the 
comparison of many theoretical models of frost growth. No data on frost surface temperature, 
frost density and heat flux were reported. 

Figure 4a illustrates the comparison of frost growth for changing wall temperatures. 
Three sets of data are considered here: Run #13: T =295.1 K, T =245.5 K, u =2.96 m/s, 

<o m =0.0108 kg v /kg a ( <f> v = 65.8%); Run #14: T a =294.9 K, T w = 254.9 K, u„ = 2.57 m/s, 

(o n =0.0106 kgv/kga =64.4%); Run #15 7„ =295.8 K, T n . =264.7 K, u w =2.92 m/s. 
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a)„ =0.01 17 kgv/kga =66.7%). The duration of the test is 150 min for Run #13, and 175 
min for Runs #14 and 15. Both the data and the model exhibit a monotonic increase in frost 
thickness and a decrease in frost growth rate with elapsed time. The relatively high frost 
growth rate is attributed to a thinner frost layer endowed with a lower thermal resistance. The 
model describes the experimental trend that the frost thickness increases as the wall 
temperature decreases as a consequence of increased heat transfer rate. The data suggest that 
at 150 min, the frost thickness is nearly doubled (1 10%) as the wall temperature increases by 
20 °K from 246 K to 265 K , as remarked by Yonko and Sepsy [4], For T w = 255 K, excellent 
agreement is achieved between the theory and the data for the entire duration. However there 
is some overprediction at T w = 265 K and underprediction at T w = 245 K. It is believed that 

uncertainties in the frost density model in the crystal growth period (where the frost density 
is low) could lead to errors in the predicted frost thickness, even though the growth rates in 
the DBG (Densification and Bulk Growth) period are satisfactorily predicted. The predicted 
frost surface temperature and frost density are found to increase with an increase in wall 
temperature, as is to be expected. 

Figure 4b shows a comparison of the model with measurements to reflect the effect of 
relative humidity. Two sets of data are considered here: Run #7: T m =294 K, T w = 245 K, 

= 1 .26 m/s, co x =0.009 kg v /kg a = 55.2%); Run #10: T a = 296 K, T w = 245 K, w M = 
1.37 m/s, 0) a =0.0141 kg v /kg a = 78.1%). The test duration for Run #7 is 200 min, and 
that for Run #10 is 140 min. The model is shown to describe the experimental trend of 
increasing frost thickness with increasing humidity ratio as a result of increased mass transfer 
rate. Satisfactory agreement between the prediction and the data is noticed for the two 
humidity ratios considered. As remarked in Yonko and Sepsy [4], a 42 % increase in the 
relative humidity (from 55.2 to 78.1%) increased the frost thickness by roughly 28%. The 
predicted surface temperature and frost density are found to increase with an increase in the 
humidity ratio, as is to be expected. 

A comparison of the model predictions with the data for frost thickness is presented in 
Fig. 4c to highlight the influence of air velocity. Two sets of data are considered here for 
comparison: Run #9: T„ =294 K, T w =245 K, u n = 1.28 m/s, (O m =0.0010 kg v /kg a (^ = 

65.2%); Run #15: T m =296 K, T w =265 K, =2.92 m/s, co n =0.01 17 kg v /kg a (^ = 
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66.7%). The data indicate that, for the conditions considered, air velocity on frost thickness 
has no appreciable effect on frost thickness as the air velocity is nearly doubled. The data 
also reveal that (excepting for the early crystal growth period) there appears to be a 
crossover point in time for the frost thickness as the air velocity increases: as the air velocity 
is increased, the frost thickness is initially reduced, but later slightly increased at large times. 
This behavior is clearly reproduced by the present model. Yonko and Sepsy [4] point out that 
although the data for frost thickness at the two air velocities are nearly identical, the frost 
densities for these test runs were quite different, although the density data are not reported. 
Physically, higher air velocity enhances heat and mass transfer to the frost surface. An 
increase in frost density with a nearly constant frost thickness implies that more vapor mass 
is diffused through the frost layer, thus densifying the frost layer. 

In order to clarify this matter theoretically the role of increased air velocity on frost 
densification, Fig. 4d shows the predicted frost densities for Run #9 and 15. The model 
indeed predicts that the frost densities are considerably different for the two air velocities 
under consideration, in accordance with the observations of Yonko and Sepsy [4]. To the 
author's knowledge, such a clarification with regard to the role of air velocity on frost 
densification has not been previously exemplified well from a theoretical point of view. 

From the foregoing comparisons, it is evident that the plate surface temperature has larger 
effect on frost thickness than the relative humidity, while both the plate temperature and the 
relative humidity have greater effect than the air velocity. These experimental facts are 
satisfactorily described by the present model. 

4.2. Comparison with Data of Cheng and Wit [8] 

The effects of air velocity and air temperature on frost thickness and frost density are 
illuminated in Fig. 5 by a comparison of the model predictions with the experimental data of 
Cheng and Wu [8], The square copper plate has dimensions of L = W =0.056 m, and a 
thickness h= 6 mm. The test duration is 30 min, and thus correspond to early stages of frost 
formation with a high resolution in time. Frost thickness is measured by a microscopic image 
system (registering data every five seconds). The uncertainty in the measured frost thickness 
is stated as 0.007 mm (0.3%). No data on frost surface temperature, frost density, and heat 
flux are reported. 
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Figure 5a shows the test data of Cheng and Wu [8] exhibiting the effect of air velocity 
on the frost thickness. Two sets of data are chosen here: Set #8a: T a =300.8 K, T w = 262.8 

K ( = 4.2 m/s, & =42%; Set #8b: T m = 300.8 K, T w = 262.8 K, u m = 6.0 m/s, <j> m =42%. 

The data suggest that increasing the air velocity (by 43% from 4.2 m/s to 6 m/s) does not 
have an appreciable effect on frost thickness. Cheng and Wu [8] report that photographic 
(visual) examination indicated that denser frost prevailed at increased velocity. These 
experimental trends concerning the role of air velocity on frost thickness and frost density 
confirm those embodied in the data of Yonko and Sepsy [4], see Fig. 4c. The model 
predictions shown in Fig. 5b are in accord with the test observations of Cheng and Wu [8] 
with regard to the frost thickness and frost density (model predicts at t = 30 min a density 
increase from 250 kg/m 3 to 300 kg/m 3 as the air velocity increased from 4. 2 m/s to 6.0 m/s). 

Figure 5c presents the test data of Cheng and Wu [8] exhibiting the effect of air 
temperature on the frost thickness. Two sets of data are chosen here: Set #6a: T a =296.5 K, 
T w = 261 .9 K, u a = 2.3 m/s, <p n = 76%; Set #6b: T m = 301 .6 K, T w = 261 .9 K, u„ = 2.3 m/s, 
<j> n =16%. The data suggest the frost thickness is rather insensitive to changes in air 
temperature (increase of 5 K). Cheng and Wu [8] again report that photographic (visual) 
examination revealed that denser frost prevailed at increased air temperature. The model 
predictions shown in Fig. 5d suggest a slight increase in frost thickness (a 12% increase in 
frost thickness from 1 .77 mm to 1 .98 mm) with increased air temperature (a 23% increase in 
air temperature relative to the melting temperature). The model concurs with the test 
observations with regard to the frost density , and suggests that at t = 30 min, model predicts 
a density increase from 236 kg/m 3 to 300 kg/m 3 as the air temperature is increased from 
296.5 K to 301.6 K. Physically, for a constant relative humidity of air (<j> n = constant), an 
increase in air temperature manifests itself in an increase in the moisture content of air 
(Cheng and Wu [8]). This situation will likely lead to increases in frost thickness and frost 
density, as the theory suggests. 

Figure 5e shows a comparison of the model predictions with the measurements of Cheng 
and Wu [8] for frost thickness, indicating the combined role of relative humidity and air 
velocity. Two sets of data were considered: Set #1 1 : T m = 300 K, T w = 266 1C, <j>. r = 41%, and 

- 4.2 m/s; Set #12: T a = 299.7 K, T w = 268 K, <j> r - 72%, =6.1 m/s. For both sets of 
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data, the model somewhat underpredicts the data by nearly a constant amount, independent 
of time during the frost layer growth period. As alluded to earlier, this difficulty is believed 
to arise due to the apparent inability of the model to accurately predict the relatively rapid 
frost growth rate in the early state of frost formation characterizing the crystal growth period. 
The result suggests that the frost thickness increases as the relative humidity increases, as 
was also observed in Fig. 4c. In the case of Set #12, the data exhibit a step-wise (and cyclic) 
increase in frost thickness for times in excess of about 18 min. This behavior was ascribed to 
the occurrence of melting visually observed at the frost surface (Cheng and Wu [8]). 

When the frost surface (air-frost interface) temperature approaches the freezing 
temperature, repeated cycles of melting and freezing occur (Sami and Duong [13]). This 
results in structural changes in the frost layer that tend to increase both the density and 
thermal conductivity with time without a proportional (corresponding) increase in the frost 
thickness. This phenomenon is more likely to occur in high humidity and/or high temperature 
environments than in low humidity and low temperature environments. In the latter case, the 
frost surface temperature will rise to a steady value below the melting point, such that the 
rate of frost deposition is balanced by the rate of sublimation from the frost (Sherif et al. 
[10]). Yonko and Sepsy [4] observe that in the test runs of high humidity and high plate 
temperature the frost surface temperature reaches the melting temperature within one to two 
hours after the commencement of the test. 

Figure 5f shows the predicted variation of frost surface temperature for Set #1 1 and Set 
#12. The frost surface is initially at the cold wall temperature and increases monotonically 
with time. No melting is predicted for Set #1 1 , characterized by low humidity condition. In 
the case of Set #12 (comparatively high humidity), the model computation is terminated once 
melting occurs at about 25 min. This prediction of the time for the approach of melting 
temperature is in reasonable agreement with the data for frost thickness, exhibiting the 
stepwise and cyclic variation with time commencing at about 20 min. These results suggest 
that the model is capable of predicting the frost surface temperature satisfactorily. Since the 
effect of air velocity is relatively small (as noted in Fig. 4c), it is to be expected that the 
increase in frost thickness is primarily on account of an increase in the relative humidity, and 
thus the comparison in Fig. 5e serves to reflect primarily the role of relative humidity on 
frost thickness. 
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4. 3. Comparison with Data of Lee et al. [ 6] 

Figures 6a-6d present a comparison of the model predictions with the test data of Lee et 
al. [6] exemplifying the role of air humidity and air velocity on frost thickness and frost 
surface temperature. Two sets of conditions are considered: Set #1 : T m =298.2 K, T w = 258.2 

K, = 1 .0 m/s, <(> r = 50%, 70% and 80%. Set #2: T m =298.2 K, T w = 258.2 K, ^ = 70%, 
u n = 0.5 m/s, and 1 m/s. The plate length L = 0.3 m. The test duration is 120 min. Digital 
micrometer is used for measuring the frost thickness, and frost surface temperature is 
measured with infrared radiation surface temperature detector assuming that emissivity of the 
frost surface is 0.94. The uncertainty in measured frost thickness is reported as ±0.1 mm, 
and that in measured frost surface temperature as ± 0.5 °C. No data on frost density and heat 
flux are reported. 

Figures 6a and 6b display the comparisons of the model predictions with test data (Set 
#1) for frost thickness and frost surface temperature signifying the role of relative humidity. 
Generally reasonable agreement is achieved between the model and the measurements for the 
frost thickness at the three values of <j) m = 50, 70 and 80% (Fig. 6a), although it must be 
admitted that the model somewhat underpredicts the data. In the case of <f>„= 50%, the model 
satisfactorily predicts that data up to Z = 90 min, beyond which the data show a small jump in 
frost thickness. The model satisfactorily predicts the experimental trend that the frost 
thickness increases as the relative humidity increases, as is to be expected. Fig. 6b suggests 
that the model describes the measured frost surface temperature reasonably well except for 
some overprediction in the case of <f>^ = 50% (where the data exhibit somewhat irregular 
behavior for times below about 40 min). At ^ = 80%, the measured frost surface 
temperature increases to about -4 °C, which is relatively close to the melting temperature. 

The model predicts the correct trend that the frost surface temperature increases with 
increasing relative humidity on account of increased mass transfer rate. 

Figure 6c shows a comparison of the predictions with the test data for frost thickness 
illustrating the effect of air velocity. The measurements suggest that as the air velocity 
increases the frost thickness increases, with the data for V =0.5 m/s and V = lm/s approach 
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each other for times in excess of about 1 10 minutes. The model is seen to satisfactory 
represent the observed frost thickness. According to the model, the air velocity effect on frost 
thickness is not appreciable, as already noted in Figs. 4c, 5a and 5b.. Fig. 6d shows a 
comparison of the frost surface temperature. The model compares well with data for frost 
surface temperature for u m = 0.5 m/s and 1 m/s, considering the irregular behavior in the data 
for u x =0.5 m/s. The model predicts the correct trend that the frost surface temperature (and 
thus the frost density) increases as the air velocity increases. This is to be expected because 
as the velocity increases heat and mass transfer rates increase, resulting in frost densification, 
as alluded to earlier. 

4.4. Comparison with Data of Lee et al. [7] 

Figures 7a-7d display respectively the comparison of the predictions with the 
measurements of Lee et al. [7] for frost thickness, frost surface temperature, frost density, 
and heat flux. The test conditions are: T n =288.2 K, T w =258.2 K, co x =0.00633 kg v /kg a , 

and w M =2.5 m/s. The test plate has dimensions of L = 0.3 m and W = 0. 1 5 m. The thickness 
of the frost layer is measured by a digital thermometer, and the frost surface temperature is 
measured by an infrared thermometer. The frost density is measured by weighing the frost 
mass with the aid of a chemical balance and the known (measured) frost thickness at t - 15, 
30, 60, 90, 120, 150 and 180 min. The uncertainty in the measurements are: 5.57% in x s , 
4.36% in T fs , 6.94% in p { and 3.69% in the total heat flux q tol . The temperature and 

humidity of the inlet/outlet air in the test section are measured by type-T thermocouples and 
humidity sensors respectively (Lee et al. [7]). 

From Fig. 7a, it is evident that there is considerable deviation between the prediction and 
the data for frost thickness for times in excess of about 90 min, with a maximum deviation 1 
mm at / = 120 min (25% error), even though the predicted growth rate is seen to be 
reasonably satisfactory. This error is compared to the uncertainty in the measured frost 
thickness of 6%. It is however noticed that there is considerable agreement between the 
theory and the data in the initial growth period (/ < 30 min). 

The frost surface temperature prediction is seen to be satisfactory as evident from Fig. 7b. 
At t = 1 80 min, the model prediction is within about 1 .0 °C of the test data. As pointed out in 
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Lee et al. [7], there are very few experimental data for frost surface temperature reported in 
the literature on account of the difficulty of measurement with a reasonable accuracy. 'Hie 
surface temperature increases relatively rapidly during the early stage of frost growth where 
the frost layer grows fast, but the rate of increase in frost surface temperature diminishes with 
time as the frost growth rate decreases. 

Figure 7c displays a comparison of the predicted frost density with the measurements. 
Excellent agreement is noted between the theory and experiments during the entire test 
period. The frost density increases rapidly in the early stage of frost growth (less than about 
30 min), and the rate of increase of frost density is small for times in excess of about 120 
min. As the frost density is correlated with the frost surface temperature, a good agreement 
between the theory with the data for frost surface temperature also results in similar 
agreement in frost density prediction. 

Figure 7d presents a comparison of the frost predictions with the measurements of 
sensible heat flux, latent heat flux and total heat flux (sensible and latent). The data suggest 
that the sensible heat flux is considerably larger than the latent heat flux (by a factor of about 
two), and decreases continuously during frost formation as the surface temperature of the 
frost layer increases with its growth (Lee et al. [7]). On the other hand, the latent heat flux 
remains fairly constant with time, following a small drop at the commencement of frost 
formation. The results also suggest that the total heat flux reduces from about 500 W/m 2 at 
t = 1 5 min to about 400 W/m 2 at t = 1 80 min (a 25% drop). 

4.5. Comparison with Data of Hermes et al. [9] 

Recently Hermes et al. [9] reported measurements of frost thickness, frost density and 
frost surface temperature. The measurements were performed on a square test plate with 
dimensions of L = W =0. 1 m. Two separate sets of measurements were obtained for frost 
growth and frost density. The test conditions for the frost thickness measurement are : 7 — 
289.2 K and 295.2 K, T w = 257.2 K to 269.2 K, (f> x = 50% and 80%, and u w =0.7 m/s and 
] .0 m/s. The test duration is 1 20 min. The test conditions for frost density measurements are: 
T m =289.2 K and 295.2 K, T w = 258 K, 261 K, 265 K and 269 K, <f> x = 50% and 80%, = 

0.7 m/s, and / =60 min and 120 min. The frost thickness is measured by a digital camera in 
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conjunction with image processing. The density (mass) is measured by an electronic scale. 
The surface temperature is measured by an infrared camera. The measurement uncertainties 
are: ±5//m(±2%) in frost thickness, ± 5% in frost density, and ±0.2°C in frost temperature. 

Figure 8 shows a comparison of the predictions with the test data of Hermes et al. [9] for 
frost thickness for various wall temperatures at T m =289 K, u m =0.7 m/s and^ = 80%. The 
present model favorably compares with the measurements for various wall temperatures, 
except for T w = 257 K where the error is relatively large. Both the data and the model 

indicate that the frost growth increases with a decrease in the cold plate temperature, as is to 
be expected. It is noticed that in the early frost growth period, the predicted frost thickness 
generally exceeds the test data for relatively high wall temperatures, while for the low wall 
temperature the trend is opposite. This observation again suggests that that the early growth 
period is important for accurate frost growth prediction. 

Table l shows a comparison of the predictions with the data of Hermes et al. [9] for frost 
thickness, frost surface temperature and frost density at t =60 and 120 min, for T x =289 K, 

=0.7 m/s, ^ M = 50 % and 80 % and T w = 15,-10 and -5 °C. The results suggest that with a 

decrease in wall temperature, the frost thickness increases, while the frost surface 
temperature and the frost density decrease. This effect arises due to an augmentation in the 
supercooling degree that results in a higher overall mass flux and a lower frost density 
(Hermes [9]). It is also noticed form the data that as the relative humidity increases, the frost 
thickness increases with an accompanying increase in frost surface temperature and frost 
density. These trends are correctly reproduced by the model. In general, the model 
satisfactorily describes the data for frost growth, frost surface temperature and frost density 
for t =60 min and 120 min. 

5. Discussion 

It is pertinent to enumerate major differences between the present model and the existing 
models such as those of Jones and Parker [1], Cheng and Cheng [17], and Na and Webb [2], 

Jones and Parker [1] consider that the energy crossing the frost surface by conduction 

=-[ K(T,-T,)+m,L„] (21) 
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where //;, is expressed by Eq. (5b). The present model considers Eq. (6d) instead, thus 

differing from Eq. (21) with regard to the second term on the right hand side. According to 
Cheng and Cheng [17], the consideration of Eq. (21) in the Jones and Parker model [1] could 
overestimate the latent energy released at the frost surface and leads to an undesired increase 
in frost surface temperature. Further, Jones and Parker [1] consider the thermal conductivity 
correlation due to Brian et al. [12], 

The primary difference between the present model and that of Chang and Chang [17] lies 
in the expressions for the frost density and frost thermal conductivity, as alluded to before. 
Cheng and Cheng [17] incorporate the empirical correlation for frost density due to Hayashi 
et al. [24]: 

P f = 650exp[0 .227(7y 5 -273.15) (22) 

where p f is in kg/m3, and T ^ is in degrees Kelvin, which has a limited range of 

applicability, and depends on the frost surface temperature only. Also, in their model, the 
thermal conductivity correlation due to Brian et al. [12] is considered. Referring to the 
present model, improved formulations for frost density [23] and frost thermal conductivity 
[25] , which are validated with experimental data over a broad range of frost conditions, are 
incorporated. 

A key element of the Na and Webb model [2] is that the water vapor is supersaturated at 
the frost surface, rather than saturated. In the present model, it is assumed that the water 
vapor at the frost surface is saturated, as was done in many frost models. The importance of 
the supersaturation condition has been alluded to earlier in the section on assumptions, and 
such a condition needs to be considered in improving the accuracy of the predictions of frost 
growth and densification. Another feature of the advanced model due to Na and Webb [2] is 
that the model requires the solution of two partial differential equations for the frost thickness 
and frost density (considering frost density gradient within the frost layer). The differential 
equations require an initial condition for the frost density, which generally involves a trial 
and error approach. However in the present model, only an ordinary differential equation for 
the frost thickness is solved for, and requires no initial condition for frost density. 

Although it would be useful to present detailed comparisons between the predictions of 
the present model with those of other models, such as those of [1, 2, 17], for several sets of 
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data covering a broad range of parameters, the author believes that such an undertaking calls 
for an elaborate effort, which is outside of the scope of the present investigation. The author 
hopes to accomplish such a comparison in a future investigation. 

It is expected that the proposed model can be easily extended to include the effects of 
variable parameters (wall temperature, relative humidity of air, air temperature, and air 
velocity) to deal with important practical problems. In this connection, it may be pointed out 
that Ellgas and Pfitzner [44] recently modeled the pulsed injection of cryogenic hydrogen 
into humid air (complex configuration with unsteady flow) by coupling a simple one- 
dimensional frost model to a commercial CFD code. 

6. Conclusions 

The proposed model frost growth and densification is shown to compare well with the 
test data, and is found to satisfactorily predict the frost thickness, frost surface temperature, 
frost density, and frost heat flux within the range of test conditions considered. It is revealed 
that the plate temperature has the greatest effect on frost thickness, and that both the plate 
temperature and relative humidity have considerably larger effect on frost thickness relative 
to the air velocity and air temperature. The model also exemplifies the experimental 
observation that the air velocity does not have appreciable effect on frost thickness, even 
though the frost density increases significantly as a result of increasing air velocity. The air 
temperature slightly increases the frost thickness, while it contributes to appreciable frost 
densification. The satisfactory comparisons of the model with test data over a wide range test 
of conditions suggests that the frost density and thermal conductivity models implemented in 
the proposed model seem to be well representative of the actual conditions. 
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Figure Captions 

Fig. 1 . Schematic of porous frost growth. 

Fig. 2. Correlation for the frost porosity and frost density as a function of dimensionless frost 
surface temperature and flow Reynolds number (adapted from Kandula [23]). 

Fig. 3. Comparisons of predicted thermal conductivity of frost with the correlations of Sturm 
et al. [34], Van Deusen [35], and Ostin and Anderson [5]. 

Fig. 4. Comparison of the present model with the data of Yonko and Sepsy [4] for 
frost thickness indicating the role of relative humidity, wall temperature and air velocity, 
data for (a): Run #13:7; =295 K, co n =0.0108, u m =2.96 m/s; Run #14: T m =295 K, co x = 
0.0106, =2.57 m/s; Run# 15: 7^ =296 K, co n =0.0117, u n =2.92 m/s; data for(b): Run 

#7: T„ =294 K, T w =245 K, u n = 1.26 m/s; Run #10: T„ =296 K, T w =245 K, u m = 1.37 
m/s; data for (c): Run #9: T a =294 K, T„ =245 K, <y w =0.009. 

Fig. 5. Comparison of the present model with the data of Cheng and Wu [8] for 

frost thickness showing the effects of air velocity, air temperature; and relative humidity; (a) 

and (c) are adapted from Cheng and Wu [8]; data for (a) and (b): T m =300.8 K, T w =262.8 

K, and <f> n = 42%; data for (c) and (d): T w = 26 1 .9 FC, <f) x = 76%, and u x = 2.3 m/s; data for 
(e) and (f): Set #11, T m =300.1 K, T w = 266 K, <j) n =41%, and u B =4.2 m/s; Set #12, T m = 
299.7 K, T w =268 K, <p a =72%, and =6.1 m/s. 

Fig. 6. Comparison of the present model with the data of Lee et al. [6] for 

frost thickness and frost surface temperature at varying relative humidity and air velocity. 

Fig. 7. Comparison of the present model with the data of Lee et al. [7] for 

frost thickness, frost surface temperature, frost density, and frost heat flux; data: T n =288.2 

K, T w =258.2 K, co m =0.00633 kg v /kg a , and u a =2.5 m/s. 

Fig. 8. Comparison of the present model with the data of Hermes et al. [9] for 
frost thickness at four different wall temperatures. 

Table Captions 

Table 1 . Comparison of predictions with the test data of Hermes et al. [9] for frost parameters 
(test conditions: T a =289.2 K, w M =0. 7 m/s). 
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Table 1. Comparison of predictions with the test data of Hermes et al. [9] for frost 
parameters (test conditions: T m =289.2 K, u n =0J m/s). 
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Fig. 1 . Schematic of porous frost growth. 
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Fig. 2. Correlation for the frost porosity and frost density as a function of dimensionless frost 
surface temperature and flow Reynolds number (adapted from Kandula [23]). 
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Fig. 3. Comparisons of predicted thermal conductivity of frost with the correlations 
of Sturm et al. [34], Van Deusen [35], and Ostin and Anderson [5]; adapted from Kandula 
[25]). 
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Fig. 4. Comparison of the present model with the data of Yonko and Sepsy [4] for 
frost thickness indicating the role of relative humidity, wall temperature and air velocity, 
data for (a): Run# 13: 7^ =295 K, co a =0.0108, u„ =2.96 m/s; Run # 14: T a =295 K, (o m = 
0.0106, u„ =2.57 m/s; Run# 15: 7^ =296 K, co m =0.0117, u a =2.92 m/s; data for (b): Run 
#7: T a =294 K, T w =245 K, u„ = 1.26 m/s; Run #10: T a =296 K, T w =245 K, u n = 1.37 
m/s; data for (c): Run #9: T a =294 K, T w = 245 K, co a =0.01 10. 
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Fig. 5. Comparison of the present model with the data of Cheng and Wu [8] for 

frost thickness showing the effects of air velocity, air temperature; and relative humidity; (a) 

and (c) are adapted from Cheng and Wu [8]; data for (a) and (b): T m = 300.8 K, T w = 262.8 

K, and(ZS M =42%; data for(c) and (d): T w =261.9 K, =76%, and u„ =2.3 m/s; data for 
(e) and (0: Set #11, T m =300.1 K, T w =266.0 K, <j>„ =41%, and u m =4.2 m/s; Set #12, 

T„ =299.7 K, T w =268.2 K, <f> m =72%, and u n =6.1 m/s. 
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Fig. 6. Comparison of the present model with the data of Lee et al. [6] for frost thickness and 
frost surface temperature at varying relative humidity and air velocity. 
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Fig. 7. Comparison of the present model with the data of Lee et al. [7] for 

frost thickness, frost surface temperature, frost density, and frost heat flux; data: T n =288.2 

K, T w =258.2 K, =0.00633 kg v /kga, and =2.5 m/s. 
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Fig. 8. Comparison of the present model with the data of Hermes et al. [9] for 
frost thickness at four different wall temperatures. 


